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ABSTRACT 


The prerequisite to a dynamic analysis of an oc#an 
platform is an accurate baseline measurement of the fundamental 
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of this baseline, the author uses the Runga-Kutta numerical 
integration method to determine the influence of a non- 
linear structural element on the measurement of the 
‘fundamental flexural natural frequency of the system. The 
non-linear phenomena which was investigated was the collision 
between the vertical conductor pipes and their guide collars 
on the platform. The analysis proved that for large excita- 
tion loads the response is non-linear, but that there is an 


excitation level below which the non-linearity is negligible. 
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CHAPTER 1 
Problem Definition and Model Description 

The size of an ocean platform and the variation of 
environmental excitations on it make the determination of the 
integrity of the platform a vital concern to the ocean 
engineer. A method of integrity certification which is under 
development utilizes on-location measurement of the platform 
natural frequencies to determine if a structural failure has 
occurred. This technique requires an initial,accurate 
baseline determination of the platform natural frequencies. 
These baseline measurements will be compared to values 
obtained at later dates to establish if any changes have 
occurred. However, before a shift in natural frequency can 
be related to damage, all other sources of change, such as 
live platform load variations and marine growth, must be 
evaluated. In this paper, the author addresses the impact 
of a non-linear structural dynamic element on the measurement 
of the fundamental natural frequency. 

For example purposes, the platform to be analyzed is a 
deep ocean drilling platform with a centrally located, 
unattached, conductor which extends from the platform deck 
to the ocean floor. The purpose of the analysis is to 
determine the effect of the non-linear interaction between 
the platform and the conductor on the measured natural 
frequency of the entire system in the fundamental flexural 


mode. In order to limit the problem to a manageable size, 


ee 

several constraints must be placed on the system. Principal 
among these is the use of a periodic excitation to drive 

the system. This constraint simplifies the identificatwvon 

of the platform response frequencies. Secondly, the initial 
separation between the platform and the conductor is set to 
zero. In the author's opinion, this position presents a worst 
case condition, and it is also the most likely from engineering 
considerations of the conductor installation. The excitation 
loads were varied in order to demonstrate a breakpoint 

between the linear and non-linear behavior of the structure. 
The remaining system parameters were set equal to specific 
values based on engineering judgement and the best information 
available. The greatest uncertainty in parameter estimation 
occurred in the determination of the conductor parameters. 

A more definitive investigation of these quantities should 

be performed before further system analysis is conducted. 

The platform can be modeled as a mass, spring and damper 
system which surrounds another mass, spring and damper system 
modeling the conductor. The platform mass is free to move 
along its horizontal axis, and if sufficient excitation as 
provided, collide with the conductoy mass. This collision 
causes a reduction in the platform movement while increasing 
the motion of the conductor. The figure on the next page 


Lllustrates the model. 





x1. ea) 


platform mass 
platform spring constant 
platform viscous damping coefficient 


conductor mass 

conductor spring constant 

conductor viscous damping coeffidient 

plat form displacement from equilibrium position 


conductor displacement from equilibrium position 


FIGURE 1.1 SYSTEM SCHEMATIC DIAGRAM 
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A mathematical solution for the response of the model is derived 
in the following manner: 
1. Equation of motion of the plattcrm mova 
separate ly: 
m)X; + b,X, + k,x, = F sin wt, 
or, x1 = (1/m,) (FP sin wt =eb pao (1) 
where x, is relative to the unforced equilibrium position of the platform. 
. 2. Equation of motion of the conductor movune 
separately: 
MoX5 + boX_, tuk. xX ome 
or, X. = (1/m,) (-b,x, - k,x,) (2) 


where x, is relative to the unforced equilibrium position of the conductor. 


oy Equation of motion for the combined platform 
and conductor masses: 
(m; + mo)X, + byxX, + box, + k,x, + kox. = F sin wt 
or, X1 = (1/mi+m2)(F sin wt - (bi+b2)X1 - kixi - kox2) (3) 
In equations1,2and 3, masses. damping constants and 
spring constants are denoted by m, b and k respectively. 
Subscript 1 designates the platform and subscript 2 denotes 
the conductor. Equations 1 and 2 completely describe the 
A A SOT SS the -sy¥otem-when the platform and—conductor—are—_ —__ 
in motion separately. Equation 3 describes the response 
after a collision has occurred. However, several questions 
about the system have yet to be answered. Specifically, 
what happens at the time of collision, how is the time of | 


collision determined, and how is the time of separation 


mle 
determined? The answers to these questions provide the 
remaining equations which will define the motion of the 
platform/conductor system at all times. 

To determine when a collision occurs, the gaps between 
the platform and the conductor must be continuously determined. 
When the value of the gap distance on either side of the 
conductor goes to zero, a collision has occurred. The 


following equation defines this feature mathematically. 


en — ~— 





4. Total gap distance; 
gap, + gap, = do - di (4) 
where do is the platform aperture diameter and di 1s the 


conductor outside diameter. 


54 Gap distance at any time, t: 
gap,(t) = gap,(t-1) + (x,(t) - x,(t-1)} -(x,(t) 
aeeX sitll) (Sa) 
gap,(t) = do - di - gap,(t) (5b) 
mere {t-1) defines the time increment just prior to the 
time,t, being calculated. 


At the time of collision, energy will be lost. However, 


momentum is conserved. Therefore the velocity of the system 


can be determined as shown below. 


oe 


6. Conservation of momentum: 
m)X, = (my, + m2) x 

or, x = (m,/m, + m2) % (6) 
Since both the conductor and platform are moving at the same 
velocityg the value of x is the value of x; in equation 3. 

Finally, an equation is needed to determine if the plat- 
form and conductor have separated. The contact force can be 
used to perform this determination. The diagram below 


illustrates the contact foree Salis 


eee —a- 


~_ Biaex 


F(t) K2x2 | 


me xX 


Aonnecvwvazon 


FIGURE 1.2 CONTACT FORCE COMPONENTS 


cae Contmet force: 


a Gen box + koXxXo - MX (7) 


See 

In summary, the ocean platform has been modeled and the 
mathematical equations have been written to define the motion 
of the platform system at any point in time. The next 
chapter illustrates the process of non-dimensionalizing these 
equations which will lend them to analysis on a more general 
level. A numerical integration method of solution will then 
be presented after which the results of the analysis will 


be determined. 


aA 
CHAPTER 2 


Non-dimensionalization of the System Equatiuns of Motion 


The removal of dimensions from the system equations of 
motion performs two important functions. First, the dimension- 
less parameters of the system have a greater application. 
Secondly, systematic variation of the dimensionless parameters 
will provide data which can be compared to known standards 
when the excitation is periodic. 

Eighteen parameters are contained in equations 1 through 
6 in chapter 1. The equations can be non-dimensionalized 
by proper selection of reference quantities for length > @maee 
and time. The author defines F/k,, m, and 27/9, as the 
reference quantities for length, mass and time respect lve 
The variables in equations 1 through 6 for displacement, 
velocity, acceleration, total gap distance and time are made 
dimensionless by division and differentiation using these 
reference values. The conductor parameters of mass, viscous 
damping coefficient and spring constant are non-dimensional- 
ized by dividing by the respective platform parimeters. After 
substitution of the new variables into equations 1,2,3 and 
7, the damping coefficient of the platform is divided by 
product of two times m times w,),. This ratio is defined as 
the damping factor,z,, and is dimensionless. The contact 
force variable becomes dimensionless after all the parameters 


of equation 7 have been non-dimensionnlized. 


a1 5 = 

For clarity, the contact force variable is redefined F to 
identify it as dimensionless. The remaining parameter, 
excitation frequency, is nondimensionalized by dividing by 
the platform natural frequency. The dimensionless parameters 


are summarized in Table 2.1. 


tt ee —_— Cae _ —_—— — 


Table 2.1, Summary of Non-Dimensional Parameters 


Non-dimensional | Dimensional 
Variable Ratio 


xR Ky) 
Xo/(F/k,) 

epee Tt 1)” 
k2/k, 

m,/m, 

b2/b, 

b,/2m,W, 

W/W 

GE Kau dx) 
Weer ky) dx} 
1/(F/k,) dx, 
TCR / ea) id * xs 
1/(2m/wWo1) dt 
1/(4n2 /wh1) dt? 
F/m1g 

(do - di)/(F/k,) 





The equations of motion from chapter 1 and their dimensionless 


equivalents are presented on the following pages. 


tee 


Platform motion moving separately: 
xX, = (1/m,) (F sincgt = 0) xe) (1) 
NON-DIMENSIONALIZED 


r = 47° (sin 20yit= ©) - 4 nes (1s) 


Conductor motion moving separately: 
Xp = (1/m2) (-b2k2 - K2x2) (2) 
NON-DIMENS IONALIZED 


@ = (1/y) (4701 & - 4107 pc) (2.70) 


Combined platform / conductor motion: 
X1 = (1/m,+mz2)(F sin wt -(b,+b2)x, - kyxX1-k2x2) (3) 
NON-DIMENSIONALIZED 


r= (1/1+u)(47?(sin 2my,t - rv - pc)-4m(Syr+Wo,r)) (3.1) 


Gap distances at any time,t: 
gap; (t) = gap)(t-1) + (x2(C) =e ee 

(xi (t) - xi (t-1)) (Sa) 
gap2(t) = do - di - gap,(t) (Sb) 
NON-DIMENSIONALIZED 
gapi(t) = gapi(t-1) + (c(t) - ¢(t-1))- 


(r(t) - E(t ois - (Sat) 


D.- gapi(tT) of _ Gb) 


gap2(T) 


Conservation of momentum: 
x = (m,/m, + m,) xX, (6) 


NON-DIMENSIONALIZED 


T= 17 Gl ee (6.1) 


ae 
Equation 6.1 appears inconsistent. However, it should be 
interpretted as the velocity the moment after a collision is 
equal to a constant times the value of the velocity just 
before the collision. 
wees COntact Force 
f(t) = b,x + k,x, - m,x (7) 
Non-dimensionalized 


f = Wt, + pe - (1/402) ur (Jaa) 


The equations of motion describing the platform/conductor 
system have been non-dimensionalized and are in a form which 
lends then to analysis through the use of an integrative 


method of solution. 


ater 
CHAPTERS 


Method of Solution and Error Analysis 


Since the interaction between the platform and the 
conductor are quite non-linear, analytic methods of solution 
must be bypassed for numerical means. Of these, the Runga- 
Kutta (R-K) is most fitting for this analysis because the 
method is self-starting and easily adaptable to computer 
programming. 

Classically, the R-K method expresses the difference 
between the values of y at x(n+1l) and x(n), where x(n+1) 
andx(n) are (n+1) and (n) increments of x away from a reference 
point x,, if the increments are constant as they will biomes 
this analysis. Mathematically this seateeneme becomes: 


m 
yinel)) =) ie maak. (8) 


where the a's are constant, the ye are deterr.ined by Taytior 
series expansion and 'm' is the order of the R-K method being 
used. Cursory examination of equation 8 indicates that higher 
orders of the R-K method yield a greater number of terms for 
each value of y(n+l). It is to be shown that h..gher orders also | 
yleld greater accuracy. 

The principle error in the Runga-Kutta method is the 


truncation error which is defined: 


m+] m+2 


TT = c(h) + o(h) 


where: m is the order of the R-K method 





mt Ou 
leusecwewimcrement (xX, —- x, ) 


¢c is a constant of Aaya 


mre pmtl 


The order of the term is generally smaller than the 
term and will not be considered further. The truncation error 
for any term may be defined: 

T; a nmr ¢ (number of steps), 
where (number of steps) = constant/h. 
Therefore, 
m 


T. «<h 
i 


An estimation of the magnitude of the truncation error for a 
term can be determined by comparison of the integration values 
uSing spacing increments of h and 2h, 


T = y(x.) - yS*%= ccany™ = c(hy™2™ 
m 1 1 


T = y(x,) - y$™) = c(h)", 
m 1 | 
where (2h) and (h) distinguish between integration increment 


values. Subtraction of the second equation from the first 


yields: 
ye = e(hy™(2" - 1) or, 
1 1 
yo ay (2h) 7 (2™ - 1). 
m i rt 
After actual calculation, T = 0.0....b, where b is the mth 
m 


decimal place. 
Keying on minimizing truncation error, the author selected 
a fourth order Runga-Kutta method for use in this paper. 


The numerica: approximation equations for this method are as 


ye 
follows: 


y(ne1) - y(n) = § (Kyot 2k2 #9 2K aoe 


kK, s@het (x(a Ce 

k, = h f£(x(n) + 'sh, y(n) + 3k) 
kz; = h £(x(n) + %h, y(n) + %ko) 
k, = h £(x(n) + h, y(n) + kj) 


To perform the R-K integration for the determinationse: 
the platform and conductor displacements, it is necessary to 
transform the second order differential equations of chapter 2 
into first order differential equations. This elementary 


procedure yields the equations shown below. 


Platform Displacement (moving separately) 


je RG 
5 =e 
ane: 
ie = 41-(€ sim 2a tT = r,) - 470,7Yr, (loz 
C= ec (1 aap 


Conductor Displacement (moving separately) 


C¢, = ec 
é, =, 
¢, 2 & 
¢, = (1/u) (41¥5,¢, - 4m%pc,) (2725 


Coa rc. (2.359) 


Oe 


Combined Platform/Conductor Displacement 


Le = (1/1+p) (40*(sin uy ,T = eS pc, ) 
-47(T,T> Gat (S52) 
ici le a (5.35) 


Platform/Conductor Gap 
gap, (Tt) a eve 4g) = corn) iy CM ee), 


eGo teen rE Ure) ) CoaaZ) 


Contact Force 


fevers + pe, = (1/4n*)ur, (72) 


The computer program utilized in this paper was designed 
and written by the author. A description, user's guide and a 
printout of the program with output are presented in Appendix 


A. 


Som 
CHAPTER 4 


Analysis Results 


As it was mentioned in chapter 1, the problem has 
been constrained to maintain manageability. As the first 
iteration toward understanding the response characteristics 
of the platform system, the author set the dimensionless 
ratios of the conductor to platform mass, viscous damping 
coefficient and spring constant to constant values. The 
damping factor, €,, and the total gap distance were aiise 
made equal to constants. The excitation force parameter, FOl, 
was given three discrete values for the purpose described in 
chapter 1. The ratio of excitation frequency to plat fom 
natural frequency, Y,, was given a range of seven valUcCGuiaage 
0.5 to 5.0. This range yields a summary plot of the plattomm 
displacement which is similar in appearance to a standard 
transfer function plot. The values of the platform mass 

“and natural frequency were estimated from previous platform 
analysis experience. The damping factor was also determined 

in this manner. With the values for natural frequency and 
mass, the platform spring constant was calculated. The ratios 
between conductor and platform parameters were then estimated. 
Using these percentages, the values of conductcr mass, 

viscous damping coefficient and spring constant were calculated. 


The values of the system parameters are listed in Table 4.1, 


Sy Be 
TABLE 4.1, System Parameter Values 
Parameter Value 


.009998 

03 

il 

BS 0nG 12 0, Leow . 0 ,.3..04.5.,0 
.01 

.01,0.001,0.0001 
NOveeemlO €lb=sec*) /inch 
.07 x 10° (lb-sec’)/inch 
wer Tor Clb sec) /ft 
mrecmeeelo- §6(1b-sec)/£t 
nome 10°  hby/cdn 

.08 x 107 lb/in 

.09 radians/second 

.628 radians/second 

19.5 inches 

17.2 inches 


0 
0 
0 
0 
0 
0 
Z 
Z 
3 
3 
9 
a 
2 
0 





The system was analyzed for the platform response a the 
mameec Of Y, for each excitation force. 

The output of the computer program was two plots. 
Each plot was coded on the abscissa to denote the type of 
response that was plotted and the specific values of y, and 
femeused for that computer run. A complete description of the 
coding and the individual displacement vs. tT plots for the 
twenty-one parameter combinations are presented with their 
analyses in Anpendix B. Each analysis discusses the effect 
of the excitation frequency on the response frequency, the 
dynamic equilibrium position either positive or negative of 
the static equilibrium zero position, the steady-state 
Magnitude of the response and the primary response frequen- 


cies, 


As an example, if y, and FOl were set to values of 


Bed 
0.5 and 0.001 respectively, the resulting platform dispiaeee 
ment and fourier transform are shown in Figures 4.0 and 4.01. 
The platform response has a period of 2T which is the period 
of the excitation frequency, 1 radian/second. The dynamic 
equilibrium position is slightly negative. The Fourier 
transform yields two frequencies, the excitation frequency 


and the platform natural frequency. 


FO1=001;, 0001 


—_ 
= 
iJ e 
= 
Lu 
O 
< 
3 
oO 
oP) 
a) 
= 
a 
O 
ue 
= 
< 
= 
o. 
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Figure 4.1, Platform Displacementua ea 
In Figure 4.1 the dimensionless platform displacement 
is plotted against the ratio of excitation frequency to 


platform natural frequency, ¢,. The non-linearity of the 
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system is demonstrated by the reduction in magnitude as the 
excitation force is decreased from 0.01 to 0,001. The 
decrease has also produced a frequency shift in the platform 
response to the right. The coincidental displacement plots for 
the excitation force equal to 0.001 and 0.0001 imply that 

the system is linear for these load levels. Thus, there is a 
breakpoint in the excitation force level below which the 
non-linear influence of the conductor on the platform natural 
frequency is negligible. However, the non-linearity has 
reduced the magnitude of the response. 

On the following pages, Figures 4.2 through 4.8 present 
plots of the absolute value of the Fourier transform of the 
platform displacement against the ratio of response frequency 
to the platform natural frequency which has been designated 
y*. Each figure displays the displacement magnitudes for a 
Specific value of the ratio of excitation frequency to platform 
natural frequency, Y,° The congruency of the plots for FOl 
equal to 0.001 and 0.0001 over all values of Y, reinforces 
the statements made previously concerning the non-linear 


effect of the conductor. 
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Principally, the analysis has shown that for large 
values of excitation force the system is non-linear. This 
is clearly displayed by the shift in the platform natural 
frequency. ‘if the excitation load is below an as yet undeter- 
mined level, the natural frequency of the platform is unaffect- 
ed by the non-linear influence of the conductor, However, 
the non-linearity reduces the steady-state platform magnitude. 
The analysis consistently attributes the response frequency 
of the platform to a function of excitation frequency and 
the platform natural frequency. Normally, the platform 
would respond at the excitation frequency in the steady- 
state mode. However, the analysis was conducted by numerical 
integration until steady-state was attained. Therefore, when 
the Fourier transform of the platform displacement was taken, 
the transient effects of the problem "start-up" were included 
and presented in the transform results. This transform is 
more desirable than a pure steady-state transform since 
actual measurements on the platform are taken during random 
excitationand the transient effects will not have an oppor- 
tunity to die out. 

With the dimensionless excitation force equal to 0.01, 
the dynamic equilibrium of the system shifted markedly 
negative if the excitation frequency was near the platform 
natural frequency. At frequencies above the platform natural 
frequency, the dynamic equilibrium always shifted positively, 


but the shift was nearly indistinguishable. The dynamic 
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equilibrium position shift was inconsistent when’ the exci- 
tation frequency equalled the platform natural frequency. 
The period of the maximum platform displacement in steady- 
State was coincidental with the excitation period throughout 
the analysis. No further conclusions can be drawn from the 
data. 

The Runga-Kutta numerical integration method has been 
an effective means of determining the dynamic characteristics 
of the system. It is evident that this analysis has raised 
more questions than it has answered because of the imposed 
parameter constraints. Even with these limitations, if several 
more of the parameters, i.¢e., the platform damping faceam 
and the conductor/platform ratios were varied, the system 
response characteristics could be more accurately determined. 
The system should then be analyzed using an excitation approx- 
imating the environmental conditions to which the platform will 
be subjected to determine the fundamental flexural baseline 
response. 

The computer program has some minor problems which 
could be eliminated for greater efficiency. The program 
was designed to analyze the system with one set of parameter 
values per computer run. However, it would be more economical 
to have the capability to increment those parameters chosen 
for variation over a set of values. This improvement would 
keep the program on the computer until all of the variations 


had been completed. The Runga-Kutta section of the program, 


ae 
although correct and complete, does not allow for varying 

the incrementing variable. If the difference between two 
adjacent, integrated values of a variable reaches a magnitude 
that is too large, the program may overflow the computer 
memory capacity and generate output which contains only error 
messages. The author endured this perennial set-back 

because of the input requirements of the Fourier transform 


and plotting subroutines. 
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APPENDIZ A 


4th order 


The computer program was written to utilize a 
Runga-Kutta (R-K) numerical integration approximation for 
determining the dimensionless displacement of an ocean 
platform with a centrally located conductor. The values of the 
displacement are then printed out, plotted as a function of T 
and subjected to a Fourier transform subroutine whose output 
absolute values are plotted as a function of the w determined 
in the transform. Because the equations of the program are 
non-dimensionalized, the values of w used for the Fourier 
transform plot are a function of T, not real time, t. To 
convert w(fourier transform) to an w(real time), one must 
divide the product of w(fourier transform) times w(platform 
natural frequency) by 27. The program has been designed to 
Compute the platform response for one set of input parameter 
emues per program execution. Thus, if Y, were assigned a 
range of values, the user must run the program on the computer 
once for each value of y,. 

The program is initiated by setting the arrays for the 
platform displacement and velocity, the conductor displacement 
and velocity, the gap values and time equal to zero. This 
ste» is followed by individually setting the values for the 
dimtnsionless parameters. The values of the dimensional 
Meeameters for excitation force, platform spring constant, 


platform aperture diameter and the conductor outside diameter 
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are then entered into the computer to calculate the value 
of the dimensionless total gap distance. The final initial 
conditions are input when the individual gap distances are 
set and the program begins execution. 

The analysis begins by computing the first integration 
values of platform displacement and velocity, and the 
conductor displacement and velocity using the ‘masses 
separated' R-K integration equations. Gap, is then calculated 
by subtracting the net change in platform displacement from the 
net change in the conductor displacement and adding this result 
to the old value of gap,. This new gap, is subtracted from 
the total gap distance to determine the new gap,. Gapiuyanm 
gap, are checked for a zero or negative value indicatiumame 
collision between the platform and the conductor. If no 
collision has occurred, time is incremented and the integration 
takes place again. If there has been a collision, the value 
of the conductor velocity is changed to equal a new value of 
the platform velocity which is determined by conservation of 
momentum. The contact force between the platform and the 
conductor is checked for a zero or negative value. If the 
contact force is zero or negative, the program returns to the 
"masses separated’ R-K integration. If the contact ferequaee 
greater than zero, the program executes the R-K approximation 
uSing the equations of motion for the combined masses. The 


program again checks for a zero contact force and incremermee 
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the R-K approximations for separated or combined masses as 
the situation dictates until the limit of the increment count 
is reached. 

At this time the platform displacement is transferred 
into a plotting subroutine for a pictorial display of the 
result. The platform displacement array is then transformed 
into a complex number array and transferred to the Fourier 
transform subroutine(FFT). The output of the FFT is transformed 
into an array of absolute magnitudes and returned to the 
plotting subroutine for display. Upon completion of the 


plot, the program ends. 
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USER CS (GU PIE 
NOTE: This computer program has been designed for execu 
on the INTERDATA M70/M80 computers of the M.I.71. Joint 
Mechanical/Civil Engineering Computer Facility 2xclusively. 
Alterations will be required to fulfill the conditions of 


plotting and Fourier transform subroutines at other facwimeges 


The computer program 1s virtually self-contained requir- 
ing only two cards for execution per computer run once 
the parameter value cards have been installed. These cards 
provide the nomenclature for the 'x' and 'y' axes of the 
plotting subroutine. These cards are placed into the jequm 
deck as shown below in the following order: platform displace- 
ment vs. T; absolute value of platform displacement vs. w&. 

The format of these two data cards is 2A40. 

The platform/conductor system parameters are included in 
the program main body. Each parameter may be varied individ- 
ually by changing the appropriate card value. 

The program execution is controlled by submitting jem. 
card deck as follows: // JOB 

// FOR 

PROGRAM DECK (including subrousmees 
RETURN 
END 

// XEQ2 

SEL) Eas 

*BC 2866 

DATA CARDS (TWO) 

// BND 


The following pages list the program and an example 


output. 


COMPLEX PWR(512) ,WTAB(256) 

INTEGER*2 XLAB(40) 

EXTERNAL AREA 

DIMENSION GAP1BR(512) ,GAP2BR(512) ,A(2,512), XSCL(4) 

COMMON R1(512) ,R2(S12) ,C1(S12) ,C2(512) , TAU(S12) , PI, AMU,GAMMA1, PSI, 
2RHO,R1KO ,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2KO,C2K1,C2K2,C2K3,L.DAMCO1, DELTAU, F01 

SET INITIAL VALUES FOR ARRAYS AND PARAMETERS 

DO 1 I= 1,512 


R1(1) 


il 
_~ 
Need 


0 


R2(I) 


1) 
ae 


0 


C1(I) 


il 
© 


i} 
= 
| 


C2 (1) 


GAPIBR(TI) 0.0 


GAP2BR(I) = 0.0 


A(1,1) 


0.0 


AC, tye 0.0 


1 TAU(I) 


0.0 


DAMCO1 = 0.03 


PI = 3.14159625 
GAMMA1 = 1.0 

DO = 19.5 
DELTAU = .1 

FO1 = .01 


FO = 79984, 


AK1 = 90784. 
DI = 17.2 

RHO = .01 

PSI = .0099998 
AMU = .1 

J=0 


DBAR = (DO - DI)/(FC/AK1) 


GAP1BR(1) = 0.0 
GAP2BR(1) = DBAR 
WRITE(S, 30) 


30 FORMAT(T10, 'TIME',T24,'PLATFORM DISP',T43,'PLATFORM VEL',T60, 'COND 
2UCTOR DISP',T80,'CONDUCTOR VEL',T99,'GAP1 BAR',//) 
START THE PROBLEM BY INTEGRATING WITH A 4TH ORDER RUNGA-KUTTA APPR 
N= 1 

4 DO=3' Fi= N,511 
L'#: TI; ° 
J=I[I+1 


TAU(J) 


TAU(I) + DELTAU 

A(2,J) = TAU(J) 

CALL RKO 

CALL RK1 

CALL RK2 

CALL RK3 

R1(J) = AREA (R1(I) ,DELTAU, R1KO,R1K1,R1K2,R1K3) 
Aiea sR 


R2(J) = AREA (R2(1I),DELTAU,R2KO, R2K1,R2K2,R2K3) 


700 


900 


100 


40 


1(J) 


C2(J) 
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AREA (C1(I) , DELTAU,C1KO,C1K1,C1K2,C1K3) 


AREA ( C2(I) , DELTAU,C2KO,C2K1 ,C2K2 ,C2K3) 


CHECK TEE GAP DISTANCES BETWEEN THE PLATFORM AND THE CONDUCTOR 


GAPIBR(J) 


GAP2BR(J) 


IF(GAPIBR(J) .LT. 
IF(GAP1BR(J) .EQ. 
IF(GAP2BR(J) .LT. 
IF(GAP2BR(J) .EQ. 


CONT INUE 


DBAR - GAPIBR(J) 


0.0) GO TO 600 
0.0) GO TO 100 
0.0) GO TO 800 


0.0) GO TO 100 


me J .EQ. S12) GO TO 400 


WHEN A GAP GOES TO 0.0, CONSERVE 


N=L+#+1 


WRITE(S,40) TAU(N) 


FORMAT(/,T6,'CONTACT OCCURRED AT 


IF(GAPIBR(N) .EQ. 


IF(GAP2BR(N) .EQ. 


0.0) GAP2BR(N) 


0.0) GAPIBR(N) 


R2(N) =(1./(1. +AMU)) * R2(N) 


C2(N) = R2(N) 


GAPIBR(I) +(C1(J) - C1i(I)) - (R1I(J) - R1IC(I)) 


MOMENTUM 


TIME,',2X,1PE12.4,/) 


INTEGRATE THE SYSTEM WITH COMBINED 


DO 5 I= N,511 


DBAR 


DBAR 


MASSES 


CHECK THE CONTACT FORCE BEFORE INTEGRATION TO ENSURE NO SEPARATION 


FI = 2.* PI * GAMMA1 * TAU(I) 


ACCEL IS A DUMMY VARIABLE FOR THE ACCELERATION AT ANY TIME TAU 


ACCEL = (4. * PI/(1. + AMU)) * (PI * SIN(FI)-DAMCO1 * R2(I) - PS 


21 * DAMCO1 * R2(I) - PI * RHO * C1(I)) 
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FHAT IS THE VARIABLE FOR THE CONTACT FORCE 
FHAT = FO1 * PSI * DAMCO1 * R2(I)/PI + RHO * FO1 * Ci(I) - AMU * 
2FO01 “* ACCEL/ sia PP a2) 

IF( FHAT .EQ. 0.0) GO TO 15 

IF( FHAT .LT. 0.0) GO TO 15 

L=I 

J=L+i1 

TAU(J) = TAU(I) + DELTAU 

A(2,J) = TAU(J) 

CALL RKOI1 

CALL RL11 

CALL RK21 

CALL RK31 

R1(J) = AREA (R1(I) ,DELTAU,R1KO,R1K1,R1K2,R1K3) 


A(1,J) = R1(J) 


R2(J) = AREA (R2(I),DELTAU, R2KO,R2K1,R2K2,R2K3) 


C1(J) = C1(1) +(RiG) - RiGee 


C203) R2 (J) 


GAPIBR(J) = GAP1BR(I) 


GAP2BR(J) 


-* 


GO 'TO 20 


GAP2BR(I) 


15 WRITE(S,25) TAU(I) ,R1(I) ,R2(1I) ,C1(1) ,C2(1) , GAP1BR(I) 
25 FORMAT(T6,1PE12.4,T25,1PE12.4,T43,1PE12.4,T61,1PE12.4,T81,1PE12.4, 
2T97 ,1PE12.4) 


IF(J .EQ. 512) GO TO 400 


WRITE(5,11) TAU(J) 
11 FORMAT (/ ,T6,'THE CONTACT FORCE IS ZERO;MASS SEPARATION OCCURRED 
2AT TIME,'2X,1PE12.4,/) 
N=L+i1 
GO TO 4 
20 WRITE(5,10) TAU(J),R1(J),R2(J) ,C1(J) ,C2(J) , GAPIBR(J) 
10 FORMAT(T6,1PE12.4,T25,1PE12.4,T43,1PE12.4 ,T61,1PE12.4,T81,1PE12.4, 
2T97 ,1PE12.4) 
5 CONTINUE 
400 M = 0 
GO TO 53 
54 DO 50 I = 1,512 
50 PWR(I) = CMPLX(R1(I) ,0.0) 


NBITS = 9 


NEWTB 


1 

INV = 1 

CALL FFT (PWR,NBITS, INV,WTAB,NEWTB) 

THE VARIABLE 'S' IS THE INCREMENTAL FREQUENCY OVER WHICH THE 
FOURIER TRANSFORM WILL BE PLOTTED, AND THE CONSTANT IN THE 'S' 
EQUATION IS EQUAL TO '2*PI/(512.*DELTAU) ' 

S = 0.0 

DO 51 I = 1,512 4 


A(1,1) 


CABS(PWR(I)) 


A(2,1) S 
91S =S + .6136 


53 NROW = 2 


XSCL(1) 0.0 


XSCL(2) = 0.0 


XSCL(3) 


0.0 


XSCL(4) = 0.0 


NVARS = 1 


NPTS 512 
NX = 2 
MOVE = 0 
LABEL = 4 
{SCL = 1 
FTIME = 0.0 
LOOK = 1 
READ(8,45) (XLAB(I), I = 1,40) 
45 FORMAT (40A2) 
CALL PICTR(A,NROW, XLAB, XSCL,NVARS ,NPTS,NX,MOVE, LABEL, ISCL, FTIME, 
2 LOOK) 
M=M+1 


IF(M .EQ. 1) GO TO 54 


GO TO 35 
600 GAPIBR(J) = 0.0 
GAP2BR(J) = DBAR 


C1(J) = R1(J) 
GO TO 700 


800 GAP2BR(J) 


0.0 
GAP1BR(J) = DBAR 


C1(J) = R1(J) + DBAR 
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THE RKO SUBROUTINE DETERMINES THE K1 COEFFICIENTS DURING 


COMPONENT SEPARATION 

SUBROUTINE RKO 

COMMON R1(512),R2(512) ,C1(512) ,C2(512) , TAU(512) , PI, AMU, GAMMA1, PST, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2KO,C2K1,C2K2,C2K3,L,DAMCO1,F01, DELTAU 

PI2 = PI **2 


FI = 2.*PI * GAMMAI * TAU(L) 


R1KO = R2(L) 

R2KO = (4. *PI2) * (SIN(FI) - RI(L)) - 4. * PI * DAMCO1 * R2(L) 
C1KO = C2(L) 

C2KO = -(4. * PI * PSI * DAMCO1 * C2(L)/AMU) - (PI2 * 4. * RHO * 
2C1(L) /AMU) 

RETURN 


END 
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C THE RK1 SUBROUTINE DETERMINES THE K2 COEFFICIENTS DURING 


C COMPONENT SEPARATION 
SUBROUTINE RK1 
COMMON R1(512) ,R2(512) ,C1(512) ,C2(512) , TAU(512) ,PI AMU,GAMMA1 PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 


3C2KO,C2K1,C2K2 ,C2K3,L,DAMCO1,F0O1 , DELTAU 


PI2 = PI **2 

B = DELTAU / 2. 

T = TAU(L) + B 

El = R1(L) + B * R1KO 

E2 = R2(L) + B * R2KO 

E3 = C1(L) + B * C1KO 

E4 = C2(L) + B * C2KO 

FI = 2. * PI * GAMMA1 * T 

R1K1 = E2 

R2K1 = (4. * PI2) * (SIN(FI) - E1) - 4. * PI * DAMCO1 * £2 
C1K1 = E4 

C2K1 = -(4. * PI*PSI*DAMCO1*E4/AMU) - (PI2 * 4.*RHO*E3/AMU) 
RETURN 


END 
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THE RK2 SUBROUTINE DETERMINES THE K3 COEFFICIENTS DURING 


COMPONENT SEPARATION 

SUBROUTINE RK2 

COMMON R1(512) ,R2(512) ,C1(512) ,C2(512) , TAU(S12) , PI, AMU, GAMMA1 , PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1K0,C1K1,C1K2,C1K3,R2K0,R2K1,R2K2,R2K3, 
3C2K0,C2K1,C2K2 ,CZ2K3, L, DAMCOl, DELTA Gs 

B = DELTAU/2. 

PAZ SPL, TZ 

T = TAU(L) + B 


Ed 


R1(L) + B * R1K1 
E2 = R2(L) + B * R2K1 
ES = Ci (lLjece® Clia 
E4 = C2(L) + B * C2K1 


el 


i 
nN 


. * PI * GAMMAL * T 
R1K2 = E2 
R2K2 = (4. * PI2) * (SIN(FI) - El) - 4. * PI * DAMCO1*E2 


C1K2 = E4 


C2K2 -(4. * PI * PSI * DAMCO1 *E4/AMU) - (PI2*4.*RHO*E3/AMU) 
RETURN 


END 
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THE RK3 SUBROUTINE DETERMINES THE K4 COEFFICIENTS DURING 


COMPONENT SEPARATION 
SUBROUTINE RK3 

COMMON R1(512) ,R2(512) ,C1(512) ,C2(512) , TAU(S12) ,PI,AMU,GAMMA1, PST, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 


3C2KO,C2K1 ,C2K2 ,C2K3,L,DAMCO1, DELTAU, F01 


B = DELTAU 
Bi2.= PI **2 
T = TAU(L) + B 


El = RI(L) + B * R1K2 
E2 = R2(L) + B * R2K2 
E3 = C1(L) + B * C1K2 


E4 = C2(L) + B * C2K2 


FI = 2. * PI * GAMMAL * T 

R1K3 = E2 

R2K3 = (4. * PI2) * (SIN(FI) - El) - 4.*PI*DAMCO1*E2 
C1K3 = E4 

C2K3 = -(4.*PI*PSI*DAMCO1*E4/AMU) - (PI2* 4.*RHO*E3/AMU) 
RETURN 


END 
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THE RKO1 SUBROUTINE DETERMINES THE K1 COEFFICIENTS FOR THE 
COMBINED MASSES 

SUBROUTINE RKO1 

COMMON R1(512) ,R2(512) ,C1(512) ,C2(512) , TAU(S12) , PI, AMU, GAMMA1, PST, 
2RHO, R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
' 3C2KO,C2K1,C2K2,C2K3,L,DAMCO1, DELTAU, F01 


FI = 2. * PI * GAMMAL * TAU(L) 


PI2 = PI **2 
R1KO = R2(L) 
R2KO = (4.*PI/(1.+AMU)) * (PI*SIN(FI) -DAMCO1*R2(L) -PSI*DAMCO1*R2 


2(L) - PI * RHO * C1(L)) 
C1KO = C1(L) 
RETURN 


END 
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THE RK11 SUBROUTINE DETERMINES THE K2 COEFFICIENTS FOR THE 


COMBINED MASSES 

SUBROUTINE RK11 

COMMON R1(512) ,R2(512) ,C1(512) ,C2(512) , TAU(S12) , PI, AMU, GAMMAI , PST, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1 ,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2KO ,C2K1,C2K2,C2K3,L, DAMCO1 , DELTAU, FO1 


B 


DELTAU + B 


T 


TAU(L) + B 

FI = 2, * PI * GAMMAL * T 
El = R1(L) + B * R1KO 

£2 = R2(L) + B * R2KO 


E3 = C1(L) + B * C1KO 


R1K1 = E2 


R2K 1 


(4. * PI/(1.+ AMU))*(PI*SIN(FI) -DAMCO1*E2/PSI* DAMCO1*E2 
2-PI * RHO * E3) 

C1K1 = E3 

RETURN 


END 
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THE RK21 SUBROUTINE DETERMINES THE K3 COEFFICIENTS FOR THE 


COMBINED MASSES 

SUBROUTINE RK21 

COMMON R1(512) ,R2(512) ,C1(512) ,C2(512) , TAU(S12) , PI, AMU, GAMMA1 , PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R1K1,R2K2,R2K3, 


5C2ZKO ,C2K1 ,C2K2 ,C2K3,. 1, DAMCO] DEE AUL LOT 


B = DELTAU / 2. 
T = TAU(L) + B 
FI = 2. * PI * GAMMAI * T 


El = R1(L) + B * RIKI 

E2 = R2(L) +B * R2K1 

E3 = C1(L) + B * C1K1 

RIK2 = E2 

R2K2 =(4.*PI/(1.+AMU)) * (PI*SIN(FI) -DAMCO1*E2-PSI*DAMCO1*E2-P1*E1- 
2PI* RHO * E3) 

C1K2 = E3 

RETURN 


END 
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THE RK31 SUBROUTINE DETERMINES THE K4 COEFFICIENTS FOR THE 


COMBINED MASSES 

SUBROUTINE RK31 

COMMONR1 (512) ,R2 (512) ,C1(512) ,C2(512) , TAU(S12) , PI, AMU, GAMMA1, PST, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2KO ,C2K1 ,C2K2 ,C2K3,L,DAMCO1 , DELTAU, F01 

B = DELTAU 


T 


TAU(L) + B 

FI = 2. * PI * GAMMAL * T 

El = RI(L) + B * R1K2 

E2 = R2(L) + B * R2K2 

Bor= Cut) +B * CiK2 

R1K3 = E2 

R2K3 =(4.*PI/(1.+AMU)) *(PI*SIN(FI) -DAMCO1*E2-PSI*DAMCO1 *E2-PI* 
2RHO * E3) 

C1K3 = E3 

RETURN 


END 
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THE AREA FUNCTION COMPUTES: Y(N+1) - Y(N) FOR R1,R2,C1,C2 


FUNCTION AREA(B1,B2,B3,B4,BS,B6) 
AREA = Bl + B2/6. * (B3 + 2.*B4 + 2.*BS + B6) 
RETURN 


END 


Root 


CU GALS aC 


00 40000°0 
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00 40000°0 
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00 40000 °0 


CO-d00L6  9- 
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Uva LdVO 
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APPENDIX B 


In this appendix the individual platform system responses 
are presented. Each plot is serialized below the abscissa 
in a statement, * Run Code XYZ." The code is interpretted 
as follows: X = 1; the plot is the result of a R-K integra- 


eon . 


we 
" 


1 through 7; denotes the value of y, used in 
the R-K integration. 

Z= 1 through 3; denotes the value of FO1 used 

during the R-K integration. 


A complete description of the codes is listed in Table B.1. 


TABLE B.1, Computer Run Code Summary 





B-2 
Each plot is followed by an analysis of the platform respomee 
with respect to steady-state magnitude, the effect of the 
excitation frequency, the position of the dynamic equilib 
positive or negative, relative to the static equilibrium zero 
position and the primary response frequencies as determined by 
a Fourier transform of the platform response. The plots which 
were computed for FO1 equal to 0.0001 are included without 
an explanation of the plots. These plots coincided with the 
platform response plots for FO1l equal to 0.001 for all values 
of y,. Therefore no additional information could be obtained 


from them, 
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FIGURE B.1, Platform Response for FO1 = 0.01; y,= 0.5 

The steady-state platform displacement is 3.5 units. 
The period between maximum displacements is 2T which corres- 
ponds to the excitation frequency of 1 radian/second. The 
dynamic equilibrium position is slightly positive. The Fourier 
transform of the platform response yields primary frequencies 
of 0.42 and (0.84 radians/second. These frequencies correspond 
to the difference between the excitation frequency and the 


conductor natural frequency, and the first harmonic of this 


ime quency. 
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FIGURE B.2, Platform Response for FOi= 0.01; )))=siGmee 











The steady-state magnitude of the platform response is 
8.87 units. The period of the response is 1.25T corresponding 
to the excitation frequency. The dynamic equilibrium position 
1s 3.84 units negative. The primary frequencies Of theme 
sponse are the excitation frequency, the sum of the excitation 
frequency and the platform natural frequency, and an unident- 
ified combination of frequencies equal to 5.7 radians/second. 
More comprehensive data is required to determine the e€xacr 


components of this frequency. 
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FIGURE B.3, Platform Response for FO1= 0.01; y,= 1.0 
The average steady-state response magnitude is 5.19 units. 

The frequency of the displacement is equal to the excitation 
frequency. As the overall displacement moves in the negative 
direction, the individual cycle amplitudes increase over a 
four cycle period such that if the equilibrium position were 
normalized, a beat frequency would be prevalent. If this 
conclusion is true, the period over the four cycles should be 
equal to one-half the beat frequency period. or 2m/(w, - w.). 


Substituting the natural frequencies of the platform and 


B-6 
the conductor for w, and w, respectively, the period 1S )4gagum 
By inspection of the displacement plot, the period of the 
four cycles is approximately 4T. Thus, a beat trequency 
caused by the interaction of the platform and the conductor 
is a plausible explanation of the platform displacment plot. 
The dynamic equilibrium position is 2.36 units negative. The 
Fourier transform of the response indicates a primary 
frequency equal to the excitation frequency and the plagiggaa 


natural frequency. 
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FIGURE B.4, Platform Response for FOl= 0.01; y,= 1.5 
The steady-state response magnitude is 0.67 units. The 

period over one response cycle is 0.66T corresponding to the 
excitation frequency. However, the period between two adjacent 
peaks of equal amplitude is 2tT. The dynamic equilibrium 
position is slightly positive. From the fourier transform of 
the response, the principle frequencies are the excitation 
frequency, the platform natural frequency and a frequency equal 
to the difference between these two frequencies. The period 


of the third frequency is 2t corresponding to the predominant 
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period between displacement peaks. 
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FIGURE B.5, Platform Response for FOi= 0,01; y,= 2.0 
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The steady-state response magnitude is 0.50 units. The 
platform oscillates alternately with a period which equates to 
the sum of the platform and conductor natural frequencies and 
a period which is twice the summation frequency. The variation 
in response period implies that the platform and conductor are 
alternately in and out-of-phase. The dynamic is slightly 
Geemerve., tithe Fourier transform indicates primary frequencies 
equal to the excitation frequency and the platform natural 


frequency. 
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FIGURE B.6, Pilacrormeheenonce for FOl= 0.01; Y.= 3.0 
The average steady-state magnitude is 0.2 units. In 
steady-state, the platform response has a beat frequency 
Superimposed. The beat period is 20tT which is a frequency 
of 0.1 radians/second. The dynamic equilibrium positions. 
Slightly positive. The primary response frequencies are the 


excitation frequency and the platform natural frequency. 
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PeouRE 6.7, Platform Response for FOl= 0.01; y,= 5.0 
The steady-state magnitude is 0.012 units. Uniquely, 

the system settled into an unstable equilibrium position 
before reaching its final, steady-state equilibrium position. 
The system cycles at a frequency of 2.65 radians/second 
which equates approximately to the sum of the platform and 
the conductor natural frequencies. It appears that at this 
high excitation frequency the relative velocities of the 
conductor are at a level such that the platform and conductor 


momenta are nearly equal. When the conductor and platform 


B-12 
are out-of-phase, the momentum of the conductor is sufficient 
enough to cancel the momentum of the platform causing a 
reduction in platform motion. When the two components are 
in phase, the platform displacement nearly doubles, The 
dynamic equilibrium position coincides with the static 
equilibrium position. The Fourier transform produces frequen- 
cies from the platform response equal to the conductor nz«tural 
frequency, a frequency which is twice the conductor naguaa 
frequency and a frequency equal to the sum of the platforn 


and the conductor natural frequencies. 
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PIGURE B.8, Piatform Response for FOl= 0.001; y,= 0.5 
The steady-state magnitude is 1.35 Sots. i iempe rod 
of the response is 2T which equates to the excitation 
frequency. The dynamic equilibrium position is slightly 
negative. From the Fourier transform, the primary response 
frequencies are the excitation frequency and the platform 


natural frequency. 
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FIGURE B.9, Platform Response for FOl= 0,001; )) = Gm 
The average steady-state magnitude is 2.6 units. The 

individual displacement cycles have a period of 1.257 whien 
is the excitation frequency. There is a slight beat frequency 
superimposed which has a period of 10T and a frequency ot 
0.2 radians/second. The dynamic equilibrium positon 
slightly positive. The Fourier transform yields primary 
frequencies equal to the excitation frequency and the platform 


natural frequency from the platform response. 
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FIGURE B.10, Platform Response for FOl= 0.001; Y,= 1.90 

The average steady-state magnitude is 5.3 units. The 
platform displacement period is 11 corresponding to the 
excitation frequency and the platform natural frequency. 
A beat frequency is more porminent. The beat period is 4t 
Mien is a frequency of 0.52 radians/second. The dynamic 
equilibrium position is 0.55 units negative. From the 
Fourier transform, the primary frequency is the excitation 
frequency. Nothing in the transform indicates the component 


frequencies which make up the beat frequency. 
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FIGURE B.11, Platform Response for FOI= 0.001; 7) sles 
The average steady-state magnitude is 1.1 units. “he 


period of a displacement cycle is 0.671 which equates to the 


excitation frequency. The beat phenomena has a period of 
24t and a frequency of 0.08 radians/second. The dynamic 
equilibrium position is 0.07 units positive. The Pownues 


transform of the platform response indicates primary frequen- 
cles equal to the excitation frequency, the platform natural 
frequency and a frequency of 1.87 radians/second which reggae 


further analysis for component identification. 
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muGURE Bol2, Platform Response for FO1l= 0.001; y,= 2.0 
The steady-state displacement magnitude is 0.50 units. 

The period of the response cycle cannot be uniquely determined 
as the period varies over adjacent cycles as the displacement 
amplitude alternately increases and decreases. However, over 
four cycles the average period is 0.5tT which corresponds to the 
excitation frequency. The period of the amplitude alteration 
1s 4tT which is a frequency of 0.53 radians/second. The 

dynamic equilibrium position is negative slightly. The 


Fourier transform of the response indicates primary frequen- 
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FIGURE B.13, Platform Response for FOl= 0.001; y,= 3.0 
The average steady-state magnitude is 0.2 units. The 
displacement period is 0.331 which is the excitation frequency. 
The platform response exhibits a beat phenomena with a period 
of 24t and a frequency of 0.08 radians/second. The dynamic 
equilibrium position is slightly positive. From the Fourier 
transform of the response, the primary frequencies are the 
excitation frequency and the platform natural frequency. It 
appears that the excitation frequency is sufficiently removed 


from the platform natural frequency that excitation force 
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is an insignificant factor. This response plot « 


the response plot for FO1 equal to 0.01. 
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mouURE B.l4, Platform Response for FOl= 0.001; y,= 5.0 
This platform response plot coincides exactly with 

the plot for FO1 equal to 0.01. The steady-state magnitude 
1s 0.012 units. The dynamic equilibrium position coincides 
with the static equilibrium position. The primary frequen- 
eres from the Fourier transform of the response are the 
sum of the platform and conductor natural frequencies, the 
conductor natural frequency and a frequency which is twice 


the conductor natural frequency. 
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FIGURE B.15, Platform Response for F0Q1= 0.0001; y,= 0.5 
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FIGURE B.16, Platform Response for FOI= 0.0001; yY 0.8 


is 


fs) 
5 
4 
2 
5 
Zz 0 
= 
4 
= -2 
ee | 
oO 
ee -4 
Q 
= -6 
(ea! 
ta 
< -8 
= | 
oO 
-10 
-12 
a) 


FIGURE B.17, 
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Platform Response for FOl= 0.0001; 





ie 








0.5 : 
0.4 | 
0.3 | 


0.2 

0.4 . Hi 
0.0 | Hi 
-0.2 {! 
-0.3 










J 
1 






i 





ijl 
id 





| 


-0.4 
o s 10 1 20F 8&8 20 83 © 4 8 


ii 


PES 







Bt 





i) $ 109 6 6.65l—C iS cr! © 8&6 
TIME PARAMETER RUN CODE 173 


FIGURE B.21, Platform Response for FO1= 0.0001; y.= 
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